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ABSTRACT 

This paper considers gravitational perturbations in geometrically thin disks 
with rotation curves dominated by a central object, but with substantial contri- 
butions from magnetic pressure and tension. The treatment is general, but the 
application is to the circumstellar disks that arise during the gravitational col- 
lapse phase of star formation. We find the dispersion relation for spiral density 
waves in these generalized disks and derive the stability criterion for axisym- 
metric (m = 0) disturbances (the analog of the Toomre parameter Qt) for any 
radial distribution of the mass-to-flux ratio A. The magnetic effects work in two 
opposing directions: on one hand, magnetic tension and pressure stabilize the 
disk against gravitational collapse and fragmentation; on the other hand, they 
also lower the rotation rate making the disk more unstable. For disks around 
young stars the first effect generally dominates, so that magnetic fields allow 
disks to be stable for higher surface densities and larger total masses. These 
results indicate that magnetic fields act to suppress the formation of giant plan- 
ets through gravitational instability. Finally, even if gravitational instability can 
form a secondary body, it must lose an enormous amount of magnetic flux in 
order to become a planet; this latter requirement represents an additional con- 
straint for planet formation via gravitational instability and places a lower limit 
on the electrical resistivity. 

Subject headings: Magnetohydrodynamics(MHD) - Stars: formation- Protoplan- 
etary disks - Planets and satellites: formation 
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1. 



Introduction 



During the gravitational collapse that forms star/disk systems, magnetic fields are 
dragged in from the interstellar medium (e.g., Galli et al. 2006; Shu et al. 2006). Additional 
fields can be generated by the central star. Consideration of mean field magnetohydrody- 
namics (MHD) in these disks shows that magnetic effects produce substantial departures 
from keplerian rotation curves through both magnetic pressure and magnetic tension (Shu 
et al. 2007, hereafter S07). On the other hand, conservation of angular momentum implies 
that most of the material that eventually accretes onto the forming star initially lands on 
the disk (Cassen & Moosman 1981). As a result, disk surface densities can be high enough 
to support gravitational instability. In the limit of axisymmetric perturbations, the criterion 
for gravitational instability is determined by the value of the parameter Qt, 



where k = w~ l [d(w 2 VL) 2 / dw] l l 2 is the epicyclic frequency, Vt is the angular rotation rate, 
a is the sound speed, and £ is the surface density (Toomre 1964). In the presence of 
magnetic fields, however, the conditions required for gravitational instability are modified. 
The principal goal of this paper is to generalize the criterion of equation fll-ip to include the 
effects of magnetic fields. More specifically, we derive a generalized stability parameter Q M 
that characterizes magnetized disks. 

We note that gravitational instability can play two important roles in circumstellar disks 
during the star formation process. If the instabilities grow into the nonlinear regime, they 
can produce secondary bodies within the disk, such as brown dwarfs and giant planets. If 
the growing perturbations saturate, the gravitational torques can lead to redistribution of 
angular momentum and disk accretion. Both processes require the onset of gravitational 
instability, which is determined by the parameter Qm derived in this paper. 

The properties and evolution of magnetized disks also depend on the dimensionless 
mass-to- flux ratio A, defined by 



For example, gravitational collapse requires A > 1. As found by S07 and discussed herein, 
for realistic magnetic field strengths, this constraint inhibits the formation of giant plan- 
ets by gravitational instability in circumstellar disks. In addition, the generalized stability 
parameter Qm derived in this paper must be a function of A. 



A = 




o 



(1-2) 



As we show in this paper, the inclusion of magnetic fields leads to competing effects, some 
of which inhibit and some of which enhance gravitational instability and planet formation. 
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However, as outlined above, magnetic fields will be present within circumstellar disks. As a 
result, in order to understand disk physics, one must include the effects of magnetic fields, 
and the goal of this paper is to provide an assessment of these effects. 

This paper is organized as follows. We specify the equations of motion for magnetized 
disks in Section [2] and find their linearized counterparts in Section |3j This procedure leads 
to the dispersion relation for spiral density waves and the generalized stability parameter 
Qm- In Section @] we present numerical examples and apply the results to the observed 
protostellar source Ceph A HW2. The condition Q M > 1 is necessary for stability and implies 
a corresponding maximum disk mass, as shown in Section [5j The onset of instability and the 
derivation of Qm can be determined by setting the resistivity rj — 0; however, realistic disks 
have rj ^ and Section [H] outlines the corresponding effects of magnetic diffusion. We then 
consider giant planet formation in Section [7J In addition to deriving modified constraints on 
planet formation via gravitational instability due to magnetic effects, we find that magnetic 
disks require an additional constraint: The need to remove magnetic flux places a lower 
bound on the electrical resistivity 77. Finally, we conclude in Section [8] with a summary and 
discussion of our results. 



This section specifies the equations of motion for this problem. We include the effects 
of a poloidal magnetic field dragged into the disk during the gravitational collapse of the 
the natal cloud that produces a newly born star /disk system. This field threads vertically 
through the circumstellar disk and is pinched radially inward by viscous disk accretion. The 
accretion in these disks is believed to occur via the magneto-rotational instability (MRI; see, 
e.g., the review of Balbus & Hawley f998). In fact, an empirical formulation of the MRI 
viscosity in thin disks has been obtained by S07 using mixing length arguments. 

Consider the evolution of gas and magnetic field in a thin axisymmetric, viscously 
accreting disk of half-thickness Zq, surrounding a young star with mass M± at the origin 
of a cylindrical coordinate system (w,z). We denote the surface density of the disk by E, 
the radial velocity of accretion in the plane by u, the azimuthal velocity about the z axis 
by v, the component of the magnetic field threading vertically through the disk by B z , and 
the radial component of the magnetic field just above the disk that responds to the radial 
accretion flow by B+. The component of the Lorentz force per unit area in the plane of the 
disk can be written as 



2. Basic Equations 




(2-1) 
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where B+ = B+e m + B+e v and Vy = e^d/dm + e tp (l/m)d/d(p. The two terms in eq. ( 12- ip 
are associated to the effects of the magnetic tension and magnetic pressure (see Shu & 
Li 1997, who, however, adopt a different definition of magnetic pressure). In this paper we 
adopt the approximation 

where B z on the r.h.s. is evaluated in the midplane of the disk. In what follows, we neglect 
the toroidal component of the magnetic field threading the disk. In contrast, the stability of 
a disk with a purely azimuthal magnetic field was studied by Lynden-Bell (1996). 

With these specifications, the governing MHD equations in cylindrical coordinates in- 
clude the equation of continuity, 

<9£ 1 d , s 1 d , . . 

— + q [voYuU) + — — (Emv) = 0; 2-3 
at m dm m 2 dip 

and the radial and azimuthal components of the equation of momentum, 

E (^ + u ^ + ^_}?\_m_/dV\ + B I Bt_z^dBl 
\dt dm m dip zd J dm \dm J z=0 2tt An dm 1 



S (| + * + = _I^_I E («:) _ JL » {B » + B% (2-5) 

\ at am m a(p J m m a(p m \oip J z=0 Airm oip 

where V is the gravitational potential and II is the gas pressure integrated over the disk 
thickness. The vertical component of the induction equation takes the form 



dB z 1 
— - H 

at m 



— {mB z u) + -7r-(B z v) 
am dip 



1 d 



rjm 



B± dB, 



zq dm 



1 d ( r) dB z 



m dtp \m dip 



m dm 

' (2-6) 

where rj is the electric resistivity. Note that we ignore the viscous torque in the azimuthal 
component of the momentum equation (12-51) because the viscous timescale is much longer 
than the gravitational instability timescale. Also, the azimuthal component of the magnetic 
field in the disk that arises from the stretching of the poloidal field by differential rotation 
averages to zero in the vertical integration. 

The vacuum fields above and below the disk are treated using the Green's function 
technique (see, e.g., Shu & Li 1997). Since B is current-free outside the disk, the magnetic 
field can be derived from a scalar potential 

B = W . (2-7) 

The condition V • B = then implies that ^ satisfies Laplace's equation 

V 2 ^ = for z ^ , (2-8) 
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subject to the boundary condition 



= (2 " 9) 



The gravitational potential is given by the stellar and disk contirbutions, V = —GM^/m + Vd, 
where Vd satisfies Poisson's equation for a thin disk, 

V 2 V d = AirGE5(z). (2-10) 

For simplicity, in the following we assume an isothermal equation of state, II = a 2 S, where 
a is the local sound speed. 



3. Linearized Perturbation Equations 



This section presents a perturbation analysis of the equations of motion and derives 
a dispersion relation for spiral density waves in magnetized disks. First expand to first 
order, and use the subscript to denote zeroth order variables and 1 to denote first order 
perturbations. It is understood that every variable is evaluated in the midplane (or just 
above the midplane). We thus look for solutions that have the Fourier decomposition 

F(w,<p,t) = F {m) + Fi(c7)e i < wt - mv > , (3-1) 

where w is a complex frequency and m is a positive integer. To zeroth order, we assume the 
disk to be in a state of axisymmetric radial equilibrium: 



S dm \ dm 

where Q = vq/w. 

The first order equations become 
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(3-2) 



(3-3) 
(3-4) 

(3-5) 
(3-6) 
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The last two terms in the parenthesis in equation f)3-4p correspond to the thermal and 
magnetic pressure of the zeroth order solution, and both of these terms are much smaller 
than the magnetic tension term. As a result, we can drop them in the following analysis. To 
proceed further, we invoke the WKB approximation and write the perturbed quantities as 

F l {m) =F x {w)e ika , (3-7) 

where k is the radial wavenumber. We make the further assumption that \k\m 3> 1, i.e., 
the spiral perturbations are tightly wrapped. As a result, we can ignore all derivatives of 
the amplitude, or division by m, compared to derivatives of the phase. In other words, we 
may replace the radial derivative of a variable with multiplication by ik times its amplitude 
(for simplicity, we will omit the circumflex accent on the slowly- varying amplitudes). As an 
additional simplification, we assume l/|fc|cc7 is of the same order as a/mfl. This specification 
implies that the parameter \k\zo remains order unity. Next we note that the first order 
perturbation in surface density and gravitational potential are related by 

which follows from the asymptotic solution of the Poisson equation for Vd to the leading 
WKB order (see Shu 1992). Following a similar argument for Poisson's equation for the 
magnetic potential we also obtain 

B zl = -\k\*t, (3-9) 

where tyf is the value of \&i just above the disk (^ is an odd function of z: if B zl > 0, then 
< for z > and > for z < 0). From the definition (12- 7p we have 

dw \k\ dm \k\ 

Using the WKB approximation in conjunction with the above relationships, the first order 
equations become 

i(u -mn) }„* -iku! = 0, (3-11) 
27rGLo 

\k\Vi ( B z0 B+ Q 2 
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We analyze first the case of ideal MHD (we consider the effects of a non-zero resistivity 
in Section [6]). Setting r\ = and solving the first and last equations for u\ and B z i, we obtain 

" 1 = (a, - m!i) 2J!p (3 " 15) 

B -f s --^> < 3 - 16 > 

where the mass-to-flux ratio is defined through equation (11-21) . After substituting these 
results, the remaining equations become 

'< w - mSi)2 2Jb - 2n Vi = ika2 ^k ~ tk + ik ^ (l + l * |2o) - (3 " 17) 



■/ r^N^l / M ^ ^ m 



a 2 \k\ 2\k\z ( ^ . kB+ 
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27rGS n A 2 V A; LB 



>*0 



(3-18) 

After further algebraic manipulation, the leading order dispersion relation takes the form 

(w - mQ) 2 = k 2 - 27rGE Q |A;|e + A; 2 6a 2 , (3-19) 

where we define 



q = 1 + tt3\ and e = 1 ~^- ( 3 - 2 °) 



and e = i_I 

2vrS a 2 md 6-1 A 2 
In the limit of vanishing magnetic field, this dispersion relation reduces to the familiar form 
for spiral density waves in a gaseous disk (e.g., Shu 1992). The magnetic field threading the 
disk modifies the standard dispersion relation for an unmagnetized disk by (?) replacing the 
sound speed a by the magnetosonic speed 1 / 2 a = (a 2 + v\ ) 1 ^ 2 , where Vao = B 2 . z /2%'Eq is 
the Alfven speed at the disk midplane, and (ii) diluting the effects of gravity by a factor e if 
A > 1. If A < 1, the right hand side of equation (13- 19[) is always positive and no instability 
occurs. 

Although eq. (I3-19P is valid in genaral for a thin disk with any radial distribution of 
the mass-to-flux ratio A, it has the same form as the dispersion relation obtained by Shu 
& Li (1997) for the special case of a disk with spatially uniform A (an "isopedic" disk). In 
particular, the marginal stability of isopedic disks with A = 1 was demonstrated explicitly 
by Zweibel & Lovelace (1997). The magnetically modified Toomre Qm parameter, which 
provides the boundary of stability for axisymmetric (m = 0) perturbations, is thus given by 

(3 - 2i) 

Note that the definition of in equation f)3-20p differs slightly from that of Shu & Li (1997). 
For Q M < 1, perturbations with wavenumber between k± = k maiX (l±^/l — Q 2 M ) are unstable, 



-8 - 



with fc max = (e/Q)kj being the wavenumber of maximum growth, and where kj = 7rGSo/a 2 
is the Jeans wavenumber. Since e/0 < 1, the effect of the magnetic field is to increase the 
length scale of the gravitational instability with respect to the Jeans length scale. 

Another important factor that determines Qm in eq. f)3-2ip is the epicyclic frequency. 
Disks around young stars that have dragged in magnetic fields from the interstellar medium 
by gravitational collapse do not rotate at keplerian speeds because magnetic tension modifies 
the force balance equation (see, e.g., eq. (18) of S07 when magnetic tension dominates over 
magnetic and gas pressure). S07 showed that in magnetized disks that are viscously accreting 
by the MRI, the rotation curve is subkeplerian by a constant fraction /. In their models, the 
subkeplerian parameter / is determined by their equation (73) that states that the magnetic 
flux brought in by star formation is conserved and is left behind in the disk. Thus, for a 
given mass-to-flux ratio A, the factor / depends on the stellar mass, M* (necessary to recover 
the flux brought in by star formation), the mass accretion rate, M<j, and the system age, 
t age . For A ~ 4, S07 obtained values of / in the range 0.39-0.95 for disks around low-mass 
and massive young stars (see their Table 2). For subkeplerian disks, the epicyclic frequency 
is given by 

K = fQ K = f (GM,/zu 3 ) 1/2 . (3-22) 

Therefore, the inclusion of magnetic fields produces competing effects on the instability 
parameter Qm- The strong fields enforce subkeplerian flow, which reduces Qm and leads to 
greater instability. On the other hand, both magnetic pressure and magnetic tension act to 
increase Qm and lead to enhanced stability. 



4. Numerical Values and Observational Application 

To evaluate the numerical values of the quantities derived in the previous section, we 
write S = iituhNh, a = (3/cT/2/im^) 1//2 , k = Cl = Q, where /i is the molecular weight, Nh 
is the hydrogen column density, T is the gas temperature, and Q is the velocity gradient. 
The mass-to-flux ratio thus becomes 



It is useful to consider an observed star/disk system where these results can be applied. 
The disk around the massive protostar Cepheus A HW2 is threaded by a large scale magnetic 
field of strength B z0 « 23 mG at a radius Rq ~ 650 AU, inferred from methanol masers 
polarization (Vlemmings et al. 2010). The disk, observed in the continuum and in several 
molecular tracers (Patel et al. 2005; Jimenez-Serra et al. 2007, 2009) is seen almost edge-on, 
with the magnetic field roughly perpendicular to the disk and inclined with respect to the line 
of sight by an angle % ~ 73° (Vlemmings et al. 2010). This inclination is in agreement with the 
high aspect ratio of the CH 3 CN and SO2 integrated emission, corresponding to inclinations 
of ~ 68° and ~ 79°, respectively (Patel et al. 2005; Jimenez-Serra et al. 2007). The disk 
shows a velocity gradient Q sini « 6 km s _1 over 0.5" (Patel et al. 2005), corresponding to 
a radial range of ~ 360 AU at the distance of 725 pc. Assuming quasi-keplerian rotation, 
the velocity gradient at the radius where the magnetic field has been measured implies 
that Q sini rj 7 x 10~ 3 km s _1 AU -1 . The inferred mass of the central star then becomes 
M+ Rj 15 sin -2 i M & . The gas temperature in the outer regions of the molecular disk is 
T w 250 K (Jimenez-Serra et al. 2009), corresponding to an isothermal sound speed of a m 
1.8 /i -1 / 2 km s _1 . Assuming that the disk is thermally supported, the scale height z rj a/Vt at 
Rq is then z m 260 /i -1 / 2 sin i AU. If the density in the methanol maser region has the typical 
value nji ~ 10 9 cm -3 , the hydrogen column density is Nh ~ 7.8 x 10 24 yU -1 / 2 sini cm -2 . 

Assuming i « 73° for consistency with the magnetic field strength determination, and 
/1 = 2.33, we obtain A ~ 1.3 (in agreement with the estimate of Vlemmings et al. 2010), 
Ri 1.8, e Ri 0.45. With these values, the combined effects of magnetic pressure and tension 
increase the value of the Toomre Qt parameter by a factor of G 1 ^ 2 /e ~ 3, reinforcing the 
stability of the disk from a marginally stable Qt ~ 1.4, to a safer Qm ~ 4.3. 

Figure [T] shows the theoretical dimensionless mass-to-flux ratio A (eq. Il-2p and the 
magnetic parameter Qm (eq. 13-2 ip as a function of radius R. The disk is unstable to 
fragmentation at large radii where A > 1 and Qm < 1- The sound speed a, A, and the 
column density E (required to compute Qm) were calculated using the profiles of magnetized 
disks in S07, for a system with a central star, M* = 16 M , and mass accretion rate, 
M = 10 -4 M yr _1 , which correspond to the expected values in the circumstellar disk 
around the massive protostar protostar Cepheus A HW2. The curves correspond to different 
values of the subkeplerian constant /. The small circles indicate the values of A and Qm 
derived in the text at R = 650 AU and the best fit corresponds to / = 0.92. 

The numerical values discussed here are certainly subject to considerable uncertainties. 
We estimate that the calculated Qm could vary in the range from a minimum of 1.3 to a 
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maximum of 7.3, depending on the uncertainty on the disk's inclination, column density, 
magnetic field and velocity gradient. While it is conceivable that this rather large range will 
be reduced by future observations, this example serves as a "proof of concept" to illustrate 
the stabilizing effect of magnetic fields in circumstellar disks, and their importance for the 
process of planet formation, addressed in the following sections. 



5. Maximum Disk Mass 

In this section, we estimate the maximum disk mass that can remain stable against 
gravitational collapse in the presence of magnetic effects. We first manipulate the critical 
stability equation into the form 

where we defined the disk scale height parameter /3 = kz /cl. 

If we set Qm — 1 and consider both the flux-to- mass ratio A and /J to be a known 
functions of radius, we can solve equation (15-11) for the critical (maximum) surface density 
profile 

S max = fM , (5-2) 

where 

m = v ^{p+[P 2 + (x 2 -i) 2 } 1/2 }- ( 5 - 3 ) 

The first factor in equation (15-21) represents the critical surface density in the absence of 
magnetic effects, the second factor represents the reduction of Keplerian rotation by magnetic 
tension, and the third represents the increase in the critical mass due to magnetic support. 
Notice that M. — > 1 in the limit of unmagnetized disk (A — > oo). In general, the dimensionless 
quantities /, A, and will be functions of the radial coordinate w in the disk. However, 
one can choose constant representative values to get an idea of how disk stability depends 
on these quantities. In particular, for a thermally supported keplerian disk, the scale height 
zq ~ a/n, and thus /3 ~ 1. 

Figure [2] illustrates the manner in which the magnetic field affects the stability of a 
circumstellar disk against gravitational perturbations, as expressed by equation (15-21) . Each 
curve in the \-Qt plane, where Qt = afQx /^GY^o represents the condition Qm = 1 for 
different values of the disk scale height parameter j3. The thick curve corresponds to /3 = 1 
and the thin curves correspond to (3 = and (3 = a/2, which is the maximum value allowed 
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Fig. 1. — Values of the dimensionless mass-to-flux ratio A and magnetically modified Toomre 
Qm as function of disk radius r for the case of the circumstellar disk around the massive 
protostar Cepheus A HW2 (M* = 16 M Q , M = 10 -4 M & yr _1 ). The curves correspond to 
different values of /, where / = 0.91 (dashes), / = 0.92 (solid), and / = 0.93 (dot-dashes). 
The small circles indicate the values of A and Q M derived in the text at R — 650 AU. 
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Fig. 2. — Curves of Qm = 1 in the X-Qt parameter space for different values of the scale 
height parameter: (3 = 1 (thick curve), /3 = and /3 = V2 (th in curves). The disk is 
unstable in the region below the curves and stable above. For / = 1, the line Qt = 1 defines 
the stability boundary in the absence of magnetic effects. Note that magnetic fields act to 
increase the region for stability. 
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for a thermally supported disk (S07). For each value of j3, the portion of the plane below the 
corresponding curve represents the region of parameter space for which disks are unstable 
in the presence of magnetic effects (Qm < !)• Above the curve, Qm is larger than unity, 
and the disk is stable. The effects of magnetic pressure and magnetic tension make a disk 
more stable compared to its unmagnetized counterpart, whereas the effects of subkeplerian 
rotation destabilize the disk. For / = 1, the line Qt = 1 defines the boundary for stability 
in the absence of magnetic effects. The same curves also show the value of the inverse of the 
function M. defined in equation ( )5-3p . 

Next we define the benchmark disk mass M max integrating the critical surface density 
when there are not magnetic effects, 

M max = / d 2wdw°^ . (5-4) 

JR. nG 

As shown by many authors (starting with Adams et al. 1988), if one uses the observed 
spectral energy distributions of T Tauri star/disk systems to specify the radial distribution of 
temperature, and hence the sound speed profile a(zu), the mass scales resulting from equation 
(15-41) lie in the range M max = 0.3 - 1 M , i.e., masses comparable to those of the central stars. 
Note that this benchmark mass scale is calculated by assuming that Qt = 1 throughout the 
disk, and thus represents an upper limit on the disk mass that can be stable. In practice, 
the Toomre parameter depends on radius, Qt = Qt{^), so that much of the disk will have 
larger Qt and hence smaller surface density than used in this exercise. Considerations of 
global stability show that the maximum disk mass is lower, e.g., the maximum disk mass 
that is stable to a class of m = 1 modes is given by M d / (M* + M d ) = 3/47T (Shu et al. 1990). 
Nonetheless, this mass scale of equation (15-41) provides an interesting benchmark. For the 
case of magnetized disks, the maximum mass is increased, as shown below. 

With the inclusion of magnetic effects, the maximum disk mass that is stable to gravi- 
tational perturbations takes the form 

M max , M = J\wdw{^\ f (p^r) 2 {/3+ [/3 2 + (A 2 -l) 2 ] 1/2 } , (5-5) 

where, in general, /, A, and (3 are functions of w. In the simplest case, however, we can 
take these parameters to be constant, with representative values. In this case, the maximum 
allowed disk mass takes the form 

M maXiM = /-MM max = M max / (j^l) i 13 + + ^ ~ X ) 2 ] ■ ( 5 ~ 6 ) 

This result differs from the field- free case by a factor F = fj\4. For example, for the choice 
of parameters (/, A, (3) = (0.9, 1.3, 1), roughly corresponding to the observed values for the 
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protostellar source Cepheus A HW2 considered in the previous section, the factor is F ~ 8 
. For somewhat larger mass to flux ratios A = 2, the factor F ~ 1.8. As a result, magnetic 
fields can produce a significant increase in the allowed masses of stable circumstellar disks. 
In the limit A — > oo, we expect / — > 1 so that F — > 1. 

It is important to keep in mind that the condition Qm > 1 is necessary for the disk to 
remain stable to axisymmetric perturbations only. An important result from spiral density 
wave theory is that disks which are stable to m = perturbations can still be unstable to 
spiral perturbations with m ^ 0. As a result, equation ( I5-5P represents an upper limit to 
the maximum stable disk mass in these systems. Similarly, equation ( 15-4)) represents an 
upper limit for the maximum stable disk mass in the absence of magnetic fields. The actual 
maximum disk mass is smaller than these benchmark scales by a factor of ~ 2 (e.g., Adams 
et al. 1989, Shu et al. 1990). 



6. Inclusion of Magnetic Diffusion 

If we include the effects of magnetic diffusion, rj ^ 0, then the dispersion relation derived 
above is modified to take the form 



(u—mQ) 2 ( u — mQ — ir}- — - I = (u— mQ) 



K2 + fcV(l+ 2^)- 2 ^'( 1 -i 



- j,JM («* + k 2 a 2 _ 2 KGE \k\) - vk 2 ^^- . (6-1) 

Zq ZTlhoZQ 

A basic analysis of this expression shows that the dispersion relation has no solutions for 
which the frequency u is purely real. As a result, u must be complex and can be written in 
the form 

tu = tu + i"Y , (6-2) 

where both uq and 7 are real. In addition, since the parameter rj is expected to be small, the 
imaginary part of the frequency is expected to be much smaller than the real part, I7I uq. 
It is useful to define the functions 

A(k) = n 2 + k 2 a 2 -27rGE |fc|, (6-3) 

and 

B (k) = k 2 + k 2 a 2 ( 1 + - 2nGZ \k\ ( 1 - 1 



27rS a 2 J \ A 2 

27rGS 
A 



A(k) + —^-\k\(l + \k\z ). (6-4) 
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We also specialize to the case of axial symmetry (m = 0) so that the dispersion relation has 
the form 

u 2 (u - irj—] = cuB(k) - ii]—A(k) - V k 2 b 2 , (6-5) 

where we have also defined b 2 = B z qB^ /(2ttTiZq). Solving for the real and imaginary parts 
of the dispersion relation, and working to leading order, we find 

cu (w% + ^IV^f) = u B(k) - V b 2 k 2 , (6-6) 

and 

"o 2 f 3 ^ - = lB(k) - rj-A(k) . (6-7) 

\ z o / z o 

To leading order, Uq = B{k), so the imaginary part of the frequency (the growth rate) takes 
the simple form 

|*| 



7 = %o 



\ A{k) 



(6-8) 



B(k)_ 

Since B(k) > A(k), with equality only in the limit \k\ — > or A — > oo, the growth rate is 
always positive, so that the solutions decay like exp(— yt). 

This decay of the perturbation solutions is expected, in general terms, because the disk 
must spread for the case i] ^ 0. However, the specific form 7 ~ i]\k\/z is less obvious. 
Nonetheless, this result can be derived by solving the magnetic induction equation in the 
limit of an infinitesimally thin disk (and this calculation is carried out in the Appendix). 



7. Giant Planet Formation by Gravitational Instability 

Gravitational instabilities in circumstellar disks can, in principle, lead to the formation 
of giant planets, or somewhat larger secondary bodies such as brown dwarfs (e.g., Boss 2001). 
In disk systems with significant magnetic support, however, the formation of secondaries is 
highly suppressed, as outlined in this section. The formation of secondary bodies requires 
both the onset of gravitational instability and sufficiently short cooling time scales (e.g., 
Gammie 2001). The required compromise between these two constraints is modified by 
sufficiently strong magnetic fields and is calculated in this section (compare with Rafikov 
2005). Even when these two constraints are met, so that gravitational instability could 
in principle produce secondaries, the magnetic flux problem remains; in other words, the 
forming protoplanet must reduce it's magnetic flux in order to contract to planetary sizes. 
This latter issue is also addressed below. 
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We first consider the coupled constraints of gravitational instability and sufficient cool- 
ing. In order for gravitational perturbations to grow, the stability parameter Qm must 
be sufficiently small. Although the growth of axisymmetric instabilities requires Qm < 1, 
spiral gravitational instabilities can grow in many star/disk systems under the weaker con- 
dition Qm < q* ~ 2 (e.g., Adams et al. 1989), where the parameter g* defines the required 
threshold. This constraint, a necessary condition for gravitational instability, takes the form 

where the dimensionless parameters and e are defined above (equation |3-20j ). 

The survival of gravitational instabilities requires that the cooling time £ coo i is sufficiently 
short, and this condition takes the form flxtcooi < a CO oi, where the value of a coo i ~ 3. Note 
that this constraint uses the Keplerian rotation rate Qk', any departures from Keplerian can 
be incorporated into the value of a coo i. Following previous authors (e.g., Gammie 2001), we 
specify the form of the cooling time according to 

7-1 2aT 4 v ; 

where r is the optical depth of the disk, 7 is the adiabatic index, and T is the temperature. 
The cooling time constraint then takes the form 

where we have defined the dimensionless parameter = 2a coo i(7 — l)/( r + V r )- Note that 
for typical values a coo i = 3, 7 = 5/3, and r = 1, the parameter = 2. 

By combining the constraints of equations (17-ip and (17-3p . one can show that the disk 
temperature — at the radial location where secondary formation could occur - obeys the 
ordering 

(E) 3 f!^y > t 3 > - nxJ: (7 4) 

where we have used the relation T = fia 2 /k, and 



(-'.) 



By eliminating the temperature using the end points of equation (17-41) . one finds the following 
constraint on the surface density 

^f^y^ ^ v (7 . 6) 



kJ \ J p*a \7rG J \p, J (p*(j) l / b 
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If we specialize to the case of a solar-mass star, with hydrogen gas (/i = 2mp), and = 2, 
the surface density constraint takes the form 



For comparison, the Minimum Mass Solar Nebula (MMSN) is expected to have a much 
lower surface density at w — 1 AU, where £i ~ 4500 g cm -2 or smaller (e.g., Kuchner 
2004). If the surface density profile is a negative power-law with index p < 2.1, then 
gravitational instability can operate in the outer disk. For the parameter values listed here, 
planet formation via gravitational instability can only take place at radii exceeding the bound 



where the final equality specializes to the index p = 3/2, often used for the MMSN. Recall 
that the parameter / determines the degree to which the disk is subkeplerian and J 7 is 
defined through equation (17-51) . Although the result is a somewhat complicated function of 
disk parameters (especially those that characterize the magnetic field), the required radius is 
always large. Planet formation via gravitational instability can only take place in the outer 
regions of large disks. We note that most disks surrounding low-mass stars have outer radii 
of order 100 AU and hence will not generally extend out to 1100 AU where giant planets 
could form. In addition, any planets that form at such large distances would have trouble 
migrating inward to the locations where (most of) the current sample of extrasolar planets 
resides (a < 10 AU). This channel of planet formation, via gravitational instability, is thus 
expected to have a limited impact on current observations 

The MMSN, with its benchmark mass scale Md ~ 0.05M©, typically has an outer radius 
of only 30 AU. For disks that have the same form for their surface density as the MMSN, 
but extend out to the radial scales of equation (17-81) . the corresponding mass is larger, about 
Md ~ 0.2M Q (/J 7 ) 1 / 2 . Notice that the constraints on radial location are weaker if the index 
p for the surface density profile is smaller. However, the constraint on the mass scale is 
nearly the same. 

Next we consider the magnetic flux problem in the context of forming secondary bodies. 
Even if the disk is heavy enough to become unstable, and the cooling time is short enough 
to allow contraction, magnetic flux must be transferred out of the region where the planet 
forms. The length scale AR that contains enough mass to form a Jovian planet is given by 
the integral over a disk annulus that contains the planet mass, i.e., 




(7-7) 



zu> 67 (fF) 6/5 



10/(21-10p) 



AU « 1100 (/J 7 ) 2 AU 



(7-8) 




Y^wdw = mp, 



(7-9) 



-18 - 



where mp ~ lmj = 2 x 10 30 g. For a MMSN disk model, where £ = Y,\{vj/R\) 3 / 2 , where 
i?! = 1 AU, this length scale is given by 

where w is the radial location of the forming planet and where we have chosen the nominal 
location zu = 1100 AU as a reference radius (see equation |7-8j ). Gravitational instability 
can potentially form fragments with Jovian mass and this length scale (10 AU), which is 
much smaller than the Hill's radius, Rh = mfap/SM*) 1 ' 3 ~ 76 AU at this location. As a 
result, the fragment can remain gravitationally bound. 

Let's now consider flux freezing. The initial fragment is threaded with a magnetic field 
B , which can be written in terms of the mass to flux ratio A through the relation B « 
27rG 1//2 S/A. If we use the required values of £ from the above analysis, and for A of order 
unity, the initial field strength is about Bq rj 10~ 4 G at the large radii where gravitational 
perturbations can grow. To form a giant planet, the fragment must contract to planetary 
size scales Rp ~ 10 10 cm. Flux freezing implies that 

BqRq ~ BpR 2 p , (7-11) 

where Bp is the magnetic field strength expected on the planetary surface. For the values 
R = 10 AU and R P = few the surface field strength would be B P rs 2 x 10 4 G. As a 
result, planet formation requires flux freezing to be compromised. 

In these disks, the resistivity rj ^ plays the role of a diffusion constant for the magnetic 
field. The time required for the resistivity to remove magnetic field from a region of size I is 
thus the diffusion time t&s and is given approximately by 

tdiff ~ i 2 h ■ (7-12) 

In this case, t ~ (zqAR) 1 ^ 2 ~ 45 AU, where the disk half thickness is given by zq = Aw, 
and the aspect ratio is taken to be A = 0. 1(^/1 00 AU) 1//4 . The value of the resistivity is 
7] rj 2 x 10 20 cm 2 s _1 , required to dissipate enough magnetic flux to meet the constraints 
posed by measurements of paleomagnetism in meteorites (Shu et al. 2006). With these 
values, the diffusion time is only t diff rs 74 yr. To leading order, the magnetic field strength 
decreases exponentially with time, with timescale t^iff ■ in order to decrease the field strength 
by the required factor of 10 4 , the system needs 4 log 10 ~ 9.2 timescales, or about 680 yr. 
The expected values of the resistivity are thus high enough to allow the magnetic field to 
diffuse out of the protoplanet on a short timescale (680 yr), much shorter than the orbital 
timescale at the large radii where giant planets could form. Nonetheless, the required flux 
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loss places a lower bound on the resistivity, i.e., the diffusion time scale t^is from equation 
( I7-I2p must be shorter than the cooling time. This constraint can be written in the form 



where adis is a dimensionless parameter of order unity and the numerical estimate was 
obtained for the density of the MMSN at 1100 AU. 

To conclude this section: The loss of magnetic flux is crucial to the planet formation 
process, but this flux can be removed with reasonable values of the resistivity 77, which must 
obey the constraint of equation (I7-I3p . The bottleneck in the planet formation process is 
given by the coupled constraints of gravitational instability and short cooling timescales, 
where these constraints are modified significantly by the presence of magnetic fields. As a 
result, we expect giant planet formation via gravitational instability to be highly suppressed 
in these systems. 



This paper has generalized the dispersion relationship for spiral density waves in circum- 
stellar disks to include the leading order effects of magnetic fields for any radial distribution 
of the mass-to- flux ratio A (equation |3-19j ). This procedure results in a generalized version 
of the stability parameter (denoted here as Qm) for gravitational instabilities in magnetized 
disks, where Qm is given by equation (13-2 ip . Magnetic fields produce competing effects re- 
garding the stability of disks to gravitational perturbations: The increased pressure support 
and magnetic tension lead to greater stability; however, these same forces lead to subke- 
plerian rotation curves, which in turn lead to greater instability. The supporting terms 
generally dominate (see Figure [2]), so that magnetic effects lead to an overall suppression 
of gravitational instabilities. In particular, there exists a maximum disk mass that is sta- 
ble to gravitational perturbations. This maximum disk mass is larger for magnetized disks 
compared to those with B = (see Section \5\ and equation |5-6| ). These ideas can be 
tested through observations, as illustrated by the case of the disk surrounding the high-mass 
protostar Ceph A HW2 (Section H]). 

The inclusion of magnetic fields leads to significant modifications to the prospects for the 
formation of giant planets through gravitational instability. In order to form giant planets, 
the disk must be unstable so that Qm is small, and the cooling time must be short enough. 
Even in field-free disks, these coupled constraints limit the formation of planets to take 
place at large radii. The constraints become tighter in the presence of magnetic fields (see 




(7-13) 



8. Conclusion 



-20 - 



equations |7-6j and [7-7] ). Even in the event that gravitational instabilities do occur, and 
cooling time scales are short enough, the gravitationally bound fragments must loose large 
amounts of magnetic flux in order to contract to planetary size scales. This requirement 
leads to another constraint on the planet formation process, namely, a lower bound on the 
electrical resistivity (see equation [7-13] ) . 

The results of this paper indicate that giant planet formation via gravitational insta- 
bility is difficult and hence should occur rarely. Nonetheless, rare is not the same as never: 
Circumstellar disks that are sufficiently large (in radius) and massive could meet the con- 
straints on this paper and support secondary formation. Some of the planetary candidates 
that have discovered through direct imaging (e.g., Marois et al. 2008, Kalas et al. 2008, 
Lagrange et al. 2008) could represent examples where this process has taken place. 

This paper represents an important step toward understanding the effects of magnetic 
support in circumstellar disks. However, the generalized stability parameter Qm derived 
here applies only to the onset of gravitational instabilities, specifically, linear perturbations 
with azimuthal wavenumber m = 0. Spiral modes (with m ^ 0) can grow when axisym- 
metric perturbations are stable, so that this work should be generalized to include higher 
wavenumbers. In addition, the long term fate of gravitational perturbations depends on the 
nonlinear evolution of these magnetized disks, and hence fully time-dependent MHD calcu- 
lations should be carried out (e.g., Inutsuka et al. 2010). Since disks are expected to contain 
magnetic fields, these studies are vital to understanding both disk accretion and the possible 
formation of secondary bodies in star/disk systems. 
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A. Appendix: Decay of magnetic fields in disks 

Equation H6-8 j) shows that the exponential decay rate of a magnetic perturbation with 
radial wavenumber k in an infinitesimally thin disk is given by 

7 ~ 7] — . (Al) 
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Thus, the characteristic diffusion length scale £ in an infinitesimally thin disk is neither 
the radial wavelength of the perturbation ~ 1/k nor the disk thickness ~ z Q , but rather the 
geometric mean of the two, £ ~ (z Q /k) 1 / 2 . Since this result is curious and at variance with the 
usual assumption that £ ~ Zq in disks (e.g., see Parker 1979 and Zeldovich et al. 1983), in the 
following we consider the expression of the decay rate derived by Brauer & Radler (1988) 
and Krause (1990) for a magnetic perturbation in a disk of finite thickness (without the 
effects of self-gravity, pressure, and rotation), and then we take the limit for z /m — > 0. 

Consider a disk of finite thickness 2zo and infinite radius, and assume that t] is constant 
inside the disk. In this case, the magnetic induction equation becomes 

<9B 

— = -77V x (V x B) for \z\ < z . (A2) 
Next we assume a vacuum field outside the disk, 

V x B = for \z\ > z , (A3) 

and impose the condition of continuity of B at z = ±Zq. 

Any solenoidal field can be decomposed into a poloidal and toroidal components defined 
by the scalar functions P and T, respectively. Note that for an axisymmetric system P = 
dA v / dm, T = B v , and the terms "poloidal" and "toroidal" assume their usual meaning. 
The magnetic field thus can be written in the form 

B = V x (e 2 x VP) + e z x VT. (A4) 

In terms of the defining scalars P and T, equation (IA2I) separates into two pieces, 

OP dT 

— = r]VP and — = r]VT . (A5) 
dt dt v ; 

Equation (1A3j) thus becomes 

V 2 P = 0, and T = 0, (A6) 

The condition of continuity of B implies that P, T and dP/dz must be continuous on 
z = ±zq. The general solutions of equation ( 1A5I) and (1A6jl that are regular at the origin can 
be expressed in terms of Bessel functions of the first kind with the arguments km and k z z, 
where k and k z are the radial and vertical wavenumbers, respectively. 



P(m,ip,z,t) 



/ °°[ai(fc) sm(k z z) + &i(Jfe) cos{k z z)} J m (kw)e* m f^ t dk if \z\ < z 
ctky-^J^kwy™*-^ dk if \z\ > z , 
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and 

T(w,ip,z,t) 



j™[a 2 {k) sin(k z z) + b 2 {k) cos{k z z)} J m {kwy m ^ 1 dk if \z\ < z 
if \z\ > Zo, 



where k and k z are the radial and vertical wavenumbers, respectively, J m {kw) are Bessel 
functions of the first kind and ai, a 2 , fei, b 2 and c are functions determined by the initial 
field distribution and the conditions of continuity on z = ±z . The decay rate is 

1 = v(k 2 z +k 2 ), (A7) 

proportional to the inverse of the modulus of the vector wavenumber k = ke m + k z e. The 
condition of continuity of dP/dz at z = ±z leads to the relation 

(kl - k 2 ) tan(2k z z ) = 2k z k. (A8) 



The limit of an infinitesimally thin disk is recovered under the following ordering: zq <C 
k~ l <C k^ 1 <C w (the disk thickness is much smaller than the vertical wavelength, that in 



turn is much smaller than the radial wavelength). With these approximations, equation (IA8I) 
then gives k « k 2 z zo and equation flA7j) becomes 7 ~ rjk 2 ~ k/zo, as found in our case. 

Note that in a thin disk, the size of axially symmetric magnetic perturbations is limited 
vertically by the disk thickness, and, horizontally by the radial extent of the perturbation. 
The perturbative calculation shows that the diffusion rate is not determined by the smallest 
of these two scales (zq, the disk thickness) but by the geometric mean of the two. In a 
disk of finite size, the relevant diffusion scale is determined, as expected, by the modulus of 
the vector sum of the vertical and horizontal wavenumbers; boundary conditions impose a 
relation between the two wavenumbers that, in the thin disk limit, gives the derived result. 
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